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ABSTRACT 

The  theory  of  torsional  vibrations  of  a  circular,  hollow  cylinder 
with  a  piecewise  constant  periodic  variation  of  rigidity  modulus  and 
mass  density  is  developed  in  terms  of  Floquet  waves.  The  dispersion 
spectrum  is  shown  to  have  band  structure,  and  the  arrangement  of 
characteristic  sequence  at  the  end-points  of  the  Brillouin  zone  is 
studied.  The  problem  of  co-existence  of  periodic  solution  is  examined 


in  detail  and  the  regions  of  stability  and  lability  are  charted. 
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1.  INTRODUCTION  AND  GOVERNING  EQUATIONS 

Consider  an  infinite,  circular,  hollow  cylinder  with  inner 

radius  b  ,  outer  radius  a  and  thickness  2h  =  a-b,  a  >  b,  b  0. 

In  cylindrical  coordinates  (r,0,z),  the  domain  R  of  the  hollow, 

homogeneous  cylinder  is  defined  as  R  :  »  w^ere  R^  :  b  <  r  <  a, 

R2  :  |z|  <  «°  and  0  s  6  <  2v  .  In  the  longitudinal  z-direction  the 

domain  R2  is  the  union  of  denumerable  infinite  number  of  cells  R^ 

00 

with  identical  properties,  so  that  R„  :  U  R,  .  Each  cell  R,  con- 

L  a  a 

sists  of  the  union  of  two  layers  R^  and  R-  with  z-domains 

R^  :  0  <  z  <  £  and  R-  :  £  <  z  <  d  ,  where  d  5  (£  +  I)  is  the  lattice 
distance  of  a  unit  cell.  We  assume  that  the  material  properties  of  the 
elastic  layer  R^  are  the  modulus  of  rigidity  y  and  mass  density  p  , 
and  the  corresponding  material  properties  of  the  elastic  layer  Rj  are 
specified  by  the  parameters  y  and  p  ,  respectively.  We  thus  assume 
that  each  cell  consists  of  two  layers  of  lengths  £  and  £  ,  and 
material  properties  (y,p)  and  (y,p)  ,  with  the  additional  assumptions 
v/p  >  0  ,  y/p  >  0  and  y/y  >  0  .  We  also  assume  that  an  infinite 
number  of  these  unit  cells  are  bonded  together  at  their  common  interfaces 
| z |  *  nd  ,  n*  0,1,2,...  So  that  the  domain  R  defines  an  infinite, 
circular,  hollow,  layered  cylinder  with  periodic  structure  of  material 
constants,  with  real  period  d  . 

For  time-harmonic  torsional  waves  with  angular  frequency  u  ,  the 
displacement  components  (u^,  uz>  =0  ,  and  the  non-trivial  tangential 
component  of  displacement  u.  =  u  satisfies  in  an  open  domain  R  the 

U 

boundary  eigenvalue  problem 
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Tre  &R  =  0  ’  3R  5  ^r=a  and  r=b^  *  lzl  <  " 


where  for  a  homogeneous,  hollow  cylinder  R:b<r<a,  |z|  <  <*>  , 

ft  =  u/io  is  the  non-dimensional  frequency,  u  5  (7r/2h)c  is  the 
s  S 

lowest  thickness-shear  frequency  of  an  isotropic,  homogenous,  infinite 
plate  of  half-thickness  h  ,  and  c  =  (y/p)  ^  is  the  speed  of  the 
shear  wave  in  an  infinite,  homogeneous,  isotropic  elastic  medium  with 
shear  modulus  y  and  mass  density  p  ,  [1].* 

Using  Bernoulli's  method  of  separation  of  variables,  we  get  the 
separated  equations 


fct  s] + r)2  -  'l*  ■  °  • 

z"  +  [fh  ?]2Z  =  0  ,  |z|  <  - 


b  <  r  <  a 


z  <  00 


K2  =  (ft2-  52)  , 

where  we  assumed  the  product  solution  to  be  the  form  u(r,z)  -  R(r)Z(z), 
and  accents  indicate  differentiation  with  respect  to  coordinate  z  . 
Since 


Tre  "  a7  (u/r)  ’ 


it  follows  that  the  boundary  conditions  are 


~  (rR)  -  2R  -  0  on  3R  . 


Numbers  in  brackets  designate  References  at  end  of  paper. 
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In  addition  t  =  t  =0  and 
rr  rz 


T 


z6 


(5) 


Now  the  structured,  hollow  cylinder  consists  of  an  infinite 
number  of  unit  cells,  and  in  the  case  under  discussion,  each  cell  con¬ 
sists  of  two  layers.  Therefore,  for  the  two  layers  in  each  cell,  the 
governing  equations  are 


Z"  + 


—  £ 
2h  ^ 


z  -  0  ,  £2  =  <fi2-  K2)  ,  o  <  z  <  £ 


Z"  + 


t^h  *]  1  =  0  *  =  (tYfl)2-  k2>  .  *  <  *  <  d 


(6) 


where  £  and  5  are  the  longitudinal  wave  numbers  in  the  layers  of 
thickness  i l  and  £  ,  respectively,  and  y2  =  (y/p)(p/y).  Because 

the  cylinder  is  of  periodic  structure  with  period  d  ,  these  equations 
hold  from  cell  to  cell.  The  wave  numbers  £  and  £  have  therefore  to 
be  extended  as  periodic  functions  of  period  d  ,  where  y  is  a  piece- 
wise  constant.  Therefore,  in  principle,  we  have  a  second  order  differ¬ 
ential  equation  with  piecewise  constant  periodic  coefficient. 

The  general  theory  of  differential  equations  with  periodic  coeffi¬ 
cients  is  due  to  Floquet  [2],  and  an  account  of  the  theory  can  be  found 
in  books:  Ince  [3],  Stoker  [A],  and  Strutt  [5].  In  quantum-mechanical 
problems  use  of  Floquet's  theory  was  first  made  by  Kronig  and  Penney  [6], 
and  subsequently  reviewed  in  his  book  by  Brillouin  [7],  More  recently, 

Kaul  and  Herrmann  [8],  and  Delph  _et  £l  [9,10]  have  made  use  of  Brillouin's 
procedure  in  solving  some  of  the  important  wave  propagation  problems  in 
periodically  structured  elastic  solid. 
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However,  when  the  periodic  coefficients  are  piecewise  constant, 
Brillouin's  method  is  unnecessarily  long  and  arduous.  Meissner  has 
shown  that  when  the  coefficients  are  piecewise  constant,  an  exact  solu¬ 
tion  can  be  easily  obtained  in  terms  of  elementary  functions,  [11]. 

An  excellent  review  of  Meissner  equation  and  its  various  properties  is 
contained  in  a  recent  paper  by  Hochstadt,  [12].  The  chief  advantage  of 
the  method  is  that  the  si2e  of  the  characteristic  determinant  depends 
only  on  the  order  of  the  differential  equation  and  is  independent  of 
the  number  of  layers  in  the  unit  cell.  Thus  in  the  present  problem, 
where  we  are  dealing  with  a  second  order  system,  the  order  of  the 
characteristic  determinant  is  2x2,  irrespective  of  whether  there  are 
two  or  more  layers  in  the  cell. 

The  purpose  of  this  paper  is  therefore  twofold,  i)  to  illustrate 
use  of  Meissner's  method,  and  ii)  to  study  in  greater  detail  the 
problem  of  torsional  waves  in  a  hollow,  structured  cylinder.  In  addi¬ 
tion  to  the  discussion  of  dispersion  spectrum,  we  also  study  the  typical 
Lyapunov-Haupt  arrangement  of  the  characteristic  sequence  [13,14],  the 
problem  of  co-existence  of  periodic  solutions,  and  the  regions  of 
stability  and  lability.  Certain  critical  cases  and  properties  of  band- 
spectrum  are  not  discussed,  and  can  be  found  in  an  earlier  report  [15]. 


i 
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2.  SOLUTION  R(r)  IN  THE  RADIAL  DIRECTION 


To  find  the  radial  eigenfunctions  R^Cr)  »  we  first  solve  the 
boundary  eigenvalue  problem 

^  r  —  +  -^(S2r2-1)R  =  0  ,  8  =  (m</2h),  b<r<a,a>b 

U) 

— (rR)  -  2R  *  0  on  r  =  a  and  r  =  b  . 

This  is  obviously  a  self-adjoint,  singular  boundary  value  problem  with 
82  >,  0  and  real,  and  the  eigenfunctions  are  orthogonal  with  respect 
to  weight  r  ,  The  general  solution  is 

R(r)  -  A  |  Jx(6r)  -  j  BBY^Br)  ,  82  *  0  (2) 

where  A  and  B  are  two  arbitrary  constants.  Since  the  Wronskian 
rW(J^(6r),  Y^CBr)]  **  2/ir(^0),  it  follows  that  the  solution  is  valid 
for  all  admissible  values  of  separation  constant  B  .  In  particular, 
when  B  ■  0  ,  we  get  from  eqn.  (2) 

R(r)  -  Ar  +  B/r  ,  B  -  0  .  (3) 

The  two  coefficients  A  and  B  are  so  chosen  that  the  shear 

stress  t  .  vanishes  on  both  boundaries  3R  :  r  *  a  and  r  =  b. 
rfl 
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Using  the  general  solution  R(r)  in  the  boundary  condition  (1)2  , 
and  setting  the  determinant  of  the  coefficients  A  and  B  equal  to 
zero,  we  get  the  eigenvalue  equation 


P2  =  [J2C6a)Y2(Bb)  -  J2(Bb)Y2(Ba)]  =  0  , 


and  the  amplitude  coefficient  is  given  by 


B  :  A  *  4J2(ga)  :  irg2Y2(Boe)  . 


Since  the  boundary  value  problem  is  self-adjoint  ,  the  eigenvalue 
equation  (4)  has  simple  roots,  which  are  all  real  and  discrete  and  can 
be  arranged  in  an  increasing  order 


0  =  B2  <  8?  <  8~  < 
o  1  2 


<  6*  ,<  3*  < 

n— i  n 


Corresponding  to  every  eigenvalue  there  is  an  eigenfunction  and  these 
eigenfunctions  form  a  complete  set,  which  are  mutually  orthogonal  with 
respect  to  weight  r  .  For  every  eigenvalue  8^  we  can  find  Bn  in 
terms  of  A^  ,  and  therefore  the  radial  eigenfunction  Rn(r)  can  be 
written  in  the  form 


g-Y-  (-g-a)  [Y2(Ba)J1(Br)  -  J^&ajY^Sr) ]  ,  B2  *  0 


where  the  suffix  n  *  0,1,2,...  belonging  to  R(r),  A  and  S  has 

been  suppressed  for  brevity.  For  8  “  0  ,  the  eigenfunction  R  (r) 

o  o 


takes  the  simple  form 
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Rc(r)  =  Aor  ,  (7) 

and  represents  a  pure  torsional  mode. 

These  eigenfunctions  can  be  easily  normalized,  if  we  make  use  of 
the  algebraic  property  that  the  product  of  two  cylinder  functions 
C^(r)  and  D^(r)  can  be  written  as 

arC^D^r)  =  ■|^[r2(C1(r)D1(r)  +  Co(r)DQ(r)) 

-r(C  (r)D.(r)  +  C.(r)D  (r))]  .  (8) 

o  1  1  o 

Making  use  of  the  normalizing  condition 

a 

rR2(r)dr  »  1  , 
b 

and  the  algebraic  identity  just  mentioned,  we  find  that  corresponding 
to  the  eigenfunction  R(r)  ,  the  normalized  coefficient  A  in  eqn.  (6) 
is  given  by 

,  BY, (8a) 

(9) 

where 


C(r)  -  ~  r2[Y2(Ba){J2(8r)  -  Jo(gr) J2 (pr) }  + 

j|(Ba){Y2(pr)  -  Yo(Br)Y2(Br))  -  (10) 

J2(ea)Y2(Ba){2J1(8r)Y1(6r)  -  Jo(8r)Y2(6r)  -  Yo(br)J2 (Br) } ] . 


The  normalized  sc:  of  orthogonal  eigenfunctions  are  thus  given  by 


*"*  *  ■?  ( :~ 

R(r)  =  -7=7  |f  J1(3r)  - 


a/ Cir)  I 

Jb 


Jn  (3a) 

iV^(ea)  Yl(s,r) 


32  >,  0  .  .  (11) 


In  particular,  in  the  limiting  case  when  &  =  G  ,  we  have  a  pure 
torsional  mode 


R  (r)  = 


/a4-b' 


a  >  b 


(12) 


In  the  case  of  a  solid  cylinder  b  =  0  ,  3^  are  the  simple  roots 

of  the  eigenvalue  equation  J?(^a)  =  0  ,  and  the  normalized  set  of 
orthogonal  eigenfunctions  take  the  simple  form 


R(r) 


■j. 

a 


J1(er)/j2(ga)  , 


i2  >,  0 


RQ(r)  =  2r/a“  ,  32  =  0  . 


a  yj 


If  we  write  b/a  =  t  ,  then 


3a  =  <  = 


(1-t) 


3b 


•t; 


and  the  eigenvalue  equation  v-.;  takes 


i  o: 


*^',(^)^'i  (L,k)  —  a^(fc.c)\^(K^  u  ■ 


(14) 


For  t  <  1  ,  the  asymptotic  expansion  for  the  p-th  root  of  this  trans¬ 
cendental  equation  is  given  by  the  McManon  series  [  16]  , 


9 


Kp  :  p[l  +  30tq2(l-  3q2  (1  +  lit  +  t2) 


-  18q4(ll  -  9t  -  109t2  -  9t3  +  lit4)  +  ...}]  , 


(15) 


where 


=  U--t 

4irpt  * 


p  =  1,2,3,.., 


(16) 


Evidently,  the  lowest  root  of  eqn.  (14)  is  k  =  0  .  The  higher  roots 

o 

of  this  equation  are  given  with  sufficient  accuracy  by  the  asymptotic 
formula  (15).  For  t  =  2/3  ,  the  first  three  roots  are  =  1.030  4056, 
<2  ~  2.015  6684  ,  =  3.010  5057.  For  p  >  3  the  higher  roots  become 

increasingly  accurate.  For  t  =  1/3,  <2  *  2.115  5496,  =  3.081  0800, 

=  4.061  9010,  etc. 
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3.  DISPLACEMENT  FIELD  Z(z)  IN  ONE  CELL 

The  domain  of  a  typical  cell  R^  *  R^  is  the  union  of  two  layers 
Rf  x  R^  and  R^  x  R-  with  an  interface  R^.  .  The  displacement 

components  Z(z)  6  R^  and  Z(z)  6  R-  are  governed  by  the  ordinary 
differential  equations  (1.6) 

Z"(z)  +  a2Z(z)  =  0  , 

Z"(z)  +  a2Z(z)  =  0  , 

where 

V  r  -  _  IT  r 

a  2h  C  '  a  2h  ?  * 

£2  =  (ft2  -  k2)  ,  X2  =  ((yft)2  -  k2)  ,  (2) 

y2  =  (c/c)2  ,  c2  =  p/p  ,  c2  =  p/p  . 

The  continuity  of  displacement  u(r,z)  and  shear  »i"css  t  r 

z6 

at  the  interface  R^  nh  requires  that  the  functions  Z(.z)  and  Z(z) 
must  satisfy  the  continuity  conditions 


Z(J,)  »  Z(SL) 

pZ  '(£)  =  pZ'(i) 


(3) 


where  accent  denotes  differentiation  with  respect  to  the  coordinate  z  . 
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We  assume  the  general  solution  of  the  differential  equation 
(1)  in  the  form 


1  0 

Z(z)  =  A  cosaz  +  B  —  sinaz  ,  a^  %  0  ,  z  G 

(4) 

_  1  -0 

Z(z)  =  A  cosaz  +  B  •=■  sinaz  ,  a2  :>  0  ,  z  G  R-  . 

a  Jc 


The  two  linearly  independent  solutions 


w,  (z)  =  cosaz  ,  w.(z)  5  —  sinaz 

i  1  a 


satisfy  the  initial  conditions 


(5) 


w^o)  =  1  ,  w|(o)  =  0  ,  w2(o)  =  0  ,  w^(o)  =  1  (6) 

and  therefore  the  Wronskian  W[w^(z)  ,  w^Cz)]  =  1  .  This  establishes 
the  linear  independence  for  a2  >,  0  .  Similarly  it  follows  that  cosaz 
and  (sinaz) /a  are  also  linearly  independent  for  a2  5  0  . 

Using  the  continuity  equation  (3)  at  the  interface  z  =  l  ,  we 
can  determine  the  coefficients  A  and  B  in  terms  of  the  coefficients 
A  and  B  .  These  coefficients  are 


A  =  A(CC  +  -^  SS)  +  -(SC  -  CS)  , 
pa  a  pa 


(7) 


B 


A(CS  -  ^  SC)  +  -(SS  +  ^  CC)  , 
pa  a  pa 


a 
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where 

C  E  cosai  ,  C  S  cosai  ,  S  e  sinai  ,  S  £  sinafc  (8) 

and  the  coefficients  A  and  B  are  still  to  be  determined. 

Having  determined  the  coefficients  A  and  B  ,  the  general  solution 
in  the  fundamental  cell  U  R^  in  terms  of  the  coefficients  A  and 
B  takes  the  form 

Z (z)  =  Aw^(z)  +  Bw2<z)  ,  a2  Z  0  , 

Z(z)  =  Aw^(z)  +  Bw2(z)  ,  a2  2  0  , 

where  now 

w  (z)  E  cosaz  ,  w„(z)  E  —  sinaz  , 

A.  +•  CX 

w  (z)  e  C  cosa(z  -  Si)  -  •=§■  S  sina(z  -  i)  ,  (10) 

1  ya 

w„(z)  E  —  S  cosa(z  -  i)  +  ==■  C  sina(z  -  1)  , 

2  a  ya 

and  the  respective  Wronskians  are 

W[w1,w2]  =  1  ,  W[w1,w2J  =  y/y  .  (11) 

In  order  that  w^  and  w2  be  always  linearly  independent,  their  corres¬ 
ponding  Wronskian  must  be  non-zero.  This  requires  that  0  <  y/y  . 
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It  can  be  easily  verified  that  at  the  interface  z  =  i  of  the  two 


layers  in  a  primitive  cell 

=  wi(d)  ,  w2(i)  =  w2U) 

,  uw^(£)  =  pw^U) 


(12) 


and  therefore  the  general  solution  (9)  satisfies  the  continuity  condi¬ 
tions  (3). 

It  may  be  finally  remarked  that  this  process  can  be  easily  general¬ 
ized  if  there  are  more  than  two  layers  in  a  primitive  cell.  Such  an 
extension  to  three  and  four  layered  composite  is  contained  in  Meissner's 
paper,  [11]. 
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4.  QUASI-PERIOD IC  SOLUTIONS  IN  A  CYLINDER 
WITH  PERIODIC  STRUCTURE 


In  the  case  of  a  cylinder  with  periodic  structure  with  period  d  , 
we  require  that  at  the  common  interface  between  two  adjacent  cells,  the 
displacement  and  stress  be  continuous  and  at  least  quasi-periodic.  The 
conditions  of  continuity  at  the  cell  interface  z  *  d  are 


u(r,d  )  =  u(r,d  )  , 


^e(r*d_)  =  Tze(r’d  >  * 


Quasi-periodicity  from  cell  to  cell  requires  that 


(1) 


u(r,z+d)  =  ou(r,z)  , 
Tze(r,z+d)  =  OTz0(r,z)  , 


0  <  z  <  i 


(2) 


1 


u(r,z+d)  =  ou(r,z)  , 
^ze(r,z+d)  =  OTz0(r,z)  * 


i  <  z  <  d 


(2)2 


where  o  is  a  suitable  constant  to  be  determined,  [3],  Combining  the 
continuity  and  the  quasi-periodic  conditions  we  find  that  at  the  cell 
interface  the  appropriate  continuity  conditions  are 


u(r,d  )  -  ou(r,0+)  , 

Tze(r’d-)  “  OTze(r’0+)  ' 


(3) 
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two  roots  and  are  explicitly  given  by 
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o 


1 


o 


o2  :  l/o 


»'H?-  4] 


(13) 


We  now  distinguish  three  cases: 

i)  ] H |  <  2  .  If  | H ]  <  2,  then  from  (13)  it  follows  that  and 

are  complex  conjugate  and  of  absolute  value  unity.  In  this  case  all 
solutions  are  bounded  since  quasi-periodicity  from  cell  to  cell  requires 
that  after  n  cells 


Z(z  +  nd)  =  onZ(z),  Z(z  +  nd)  =  o°Z(z)  (14) 

where  now  |o|  =  1  . 

ii)  | H  |  >  2  .  If  |  H  j  :•  2  ,  then  from  (13)  it  follows  that  4 

and  the  two  roots  are  real.  Since  =  1  1  it  implies  that  if  one 

root  is  greater  than  unity  then  the  other  must  be  less  than  unity  in 
absolute  value.  Thus  as  the  number  of  cells  n  -*  “•  ,  one  of  the  solutions 
will  become  unbounded  and  the  other  will  decay  to  zero. 

iii)  H  =  ±  2  .  If  H  =  *  2  ,  then  =  o0  =  0  is  a  root  of  multiplicity 

2.  Since  ,,^°2  =  1  t  if  follows  that  cr  =  1  and  therefore  either  0  =  +  1 
or  o  »  -  1  ,  each  being  a  double  root.  When  0  =  1  for  H  =  2  ,  the 

solution  is  periodic  with  period  d  ,  since  in  this  case  Z(z  +  d)  =  Z(z). 

When  a  ■  -  1  for  H  =  -  2  ,  the  solution  is  periodic  with  period  2d  , 

because  in  this  case  Z(z  +  2d)  =  (-1)2  Z(z)  =  Z(z)  .  In  the  case  of 

double  roots  only  one  solution  is  basically  periodic,  the  second  linearly 
independent  solution  is  in  general  aperiodic,  [  3  ]t 

t  This  is  reminiscent  of  a  similar  situation  in  the  theory  ol  ordinary 
differential  equations,  when  the  indicial  equation  has  roots  ol  higher 
mult iplic ity . 

.  .  _  ■  -T-  *  —  — ■  — ■  —  - - 


18 


The  characteristic  equation  (7),  or  equivalently  eqn.  (11)  takes 
a  familiar  form  if  we  introduce  the  transformation 


o  =  exp2L(x  +  t)X  ,  ( i'1’  -  -  1) 


(15) 


where 


7T  Jt  -  „  It  l 

T  :  4h  *  T  ’  Ah 


(16) 


and  A  is  a  characteristic  exponent  (also  known  as  Floquet's  exponent), 
which  determines  the  phase  shift.  In  Lite  transformation  (15),  A  is 
congruent  modulo  11/(1  +  i).  For  a  given  o  t  the  value  of  A  is  deter¬ 
mined  uniquely  in  the  interval  1[X]:[Q  ,  n/(t  -t  t )].  Having  deter¬ 
mined  X  in  the  interval  I[X]  ,  all  other  values  of  X  are  congruent 
to  it  modulo  tt/(t  +  t)  . 
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5.  THE  CHARACTERISTIC  A-EQUATION 

In  terms  of  the  characteristic  exponent  X  ,  the  two  characteris¬ 
tic  equations  (4.11)2  ^  can  be  written  in  the  useful  form 

sin2 (t  +  t ) A  =  i(2  -  H)  , 

,  _  i  ( 

cos2 (x  +  t)A  =  —(2  +  H) 


The  function  H  is  defined  in  terms  of  w^(d)  and  w^(d)  which  can 
be  easily  obtained  from  w^z),  w2(z),  which  are  given  by  (3.10).  After 
some  algebraic  simplification  it  can  be  shown  that  the  two  characteristic 
equations  (1),  take  the  interesting  product  form 

sin2  (x+x) A  =  |c(tOS(tO+  S(xOC(TO  Jjc(xOS(x£)+  S(xOC(U)J  , 

(2) 

cos2  (t+t)  A  =  jc(xOC(”j-  S(tOS(U)}{c(tOC(U)-  S(xOS(t£)J  , 

where  S  and  C  are  abbreviations  for  trigonometric  sine  and  cosine 
functions,  respectively.  We  may  now  note  that  the  two  trigonometric 
functions  sin2(x+T)A  and  cos2(x+x)A  ,  are  even- functions  about  A  =  0 
and  A  «tt/[2(x+x)]  .  In  addition,  A  is  congruent  modulo  n/(x+x)  .  It 

is  therefore  sufficient  to  restrict  Floquet's  wave  number  >  to  the 
interval  I_ [ A J  :  [0,n/2(x+x)]  .  With  this  restriction  on  the  value  of 

O 

A  ,  the  frequency  spectrum  fi  versus  A  has  a  zone  structure,  and  in 
the  first  Brillouin  zone  I  ,A  varies  from  0  to  n/l2(x+x)l,  [  7  ]  . 

D 
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The  characteristic  roots  Q  of  this  characteristic  equation  can 
be  arranged  in  an  increasing  order  and  thus  can  be  indexed  by  a  suffix  n. 
It  can  be  shown  graphically  that  the  order  follows  the  well-known 
Lyapunov-Haupt  sequence,  [13,14],  In  particular,  it  can  be  exhibited 
explicitly  at  the  end-points  of  the  Brillouin  zone,  X  =  0  and 
X  =  tt/[2  (t-4-t)  ] . 

Of  particular  interest  is  the  case  when  H  =  t  2,  that  is,  when 
o  =  ±  1  .  When  o  =  1  ,  X  =  mr/(T+r)  ,  n  =  0,1,2,...  and  in  the  reduced 
zone  scheme  X  =  0  corresponds  to  the  left  end-point  of  the  interval 
I  .  Since  a  =  1  is  a  zero  of  multiplicity  2  ,  out  of  the  two  linearly 

D 

independent  solutions  only  one  is  periodic  with  period  d  ;  the  second 
solution  is  aperiodic.  Corresponding  to  the  periodic  solutions  with 
period  d  ,  the  characteristic  equation  (2)^  uncouples  into  two  equations 

1)  tan-rC  +  tanr£  =  0  , 

o*l,  H  =  2  (3) 

2)  tanrC  +  =4  tamC  =  0  . 

When  o  *  -1  ,  X  =(n+l/2)7r/(T+T) ,  n  =  0,1,2,...  and  in  the  reduced 
zone  scheme  X  =it/2(t+t)  corresponds  to  the  right  end-point  of  the 
interval  I  .  Since  a  =  -  1  is  a  zero  of  multiplicity  2,  only  one 

D 

linearly  independent  solution  is  periodic  with  period  2d.  Corresponding 
to  the  periodic  solution  with  period  2d,  the  characteristic  equation  (2)^ 
uncouples  into  two  equations 


1)  cotiC  -  ^  tam£  =  0  , 

2)  coti^  -  tant£  =  0  . 


(4) 


o  *  -1,  H  =  -2 
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One  can  easily  see  that  the  multiplicity  of  the  root  o  =  1  ,  corres¬ 
ponds  to  the  multiplicity  of  the  root  \  =  m\/(i+x)  ,  and  multiplicity 
of  the  root  o  =  -1  corresponds  to  the  multiplicity  of  the  root 
A  =(n+l/2)Tt/(T+x)  .  This  follows  from  the  fact  that  in  the  first  case 
sin2(i+i)A  and  its  first  derivative,  both  vanish  at  the  left  end-point 
of  the  Brillouin  zone.  In  the  second  case  cos2(t+t)A  and  its  first 
derivative,  both  vanish  at  the  right  end-point  of  the  Brillouin  zone. 

If  Z(z)  is  a  solution  then  Z(-z)  is  also  a  solution  because 
the  differential  equation  is  invariant  under  coordinate  reflection. 

These  two  solutions  are  in  general  linearly  independent,  unless  Z(z) 
is  an  even  or  an  odd  solution.  However,  we  have  shown  the  existence 
of  periodic  solutions  for  o  =  ±  1  .  Therefore,  either  Z(z)  is  a 
periodic  solution  or  it  is  not.  If  Z(z)  is  a  periodic  solution  then 
it  is  either  even  or  odd.  If  it  is  not  even  or  odd  then  we  can  construct 
the  solutions  [Z(z)  +  Z(-z)]  and  [Z(z)  -  Z (— z ) ]  which  are  even-periodic 
and  odd-periodic,  respectively.  In  either  case,  even-periodic  and 
odd-periodic  solutions  can  always  be  determined,  when  periodic  solutions 
exist.  Thus  at  the  left  (right)  end-point  of  the  Brillouin  zone,  even- 
and  odd-periodic  solutions  with  period  d(2d)  can  always  be  determined. 

The  dispersion  spectra  on  an  extended  zone  scheme  are  shown  in  Fig.  (1) 
for  kq  =  0  ,  and  in  Fig.  (2)  for  =  1.18920.  The  material  parameters 
are  t  -  1/3,  y  =  1/4,  p/p  =  40,  p/p  =  5/2  and  t/t  =  5.  Real  seg¬ 
ments  of  the  spectra  are  shown  as  full  lines,  imaginary  segments  are  shown 
as  dotted  lines.  In  Fig.  (2),  point  A  corresponds  to  the  limiting  case 
5  =  0  ,  and  point  B  to  the  limiting  case  when  5=0.  In  these  figures, 
hatched  regions  represent  stopping  bands. 
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6.  CHARACTERISTIC  A-EQUATION  FOR  =  0 


In  the  case  of  pure  torsional  mode  =  0  and  this  implies 
=  0,  5  =  l  .  The  characteristic  equations  (5.2)  now  take 

a  simple  form 


sin? (t+t) A  =  |c(rfi)S(rrfl)+  ~LS(m)C(yTU)  jjc(Tft) S(YTfl)+  ^S(ifl)C(yr£l)  j  , 

(1) 

cos2(t+t)A  =  |c(xfi)C(YrS2)-  •^•S(Tfi)S(YTfi)||c(Tft)C(YTSI)-  ■=^-S(Tfi)S(yTfO  . 

At  the  end-points  of  the  Brillouin  zone,  for  o  =  1  the  charac¬ 
teristic  equations  are 


1)  tanYtfi  +  tanxS!  =  0  , 

2)  tanYtfi  +  ^-=-  tanxfi  =  0  , 


o  =  1,  H  =  2 


(2) 


and  for  o  =  -1,  the  characteristic  equations  are 


1)  cotYiii  -  tanxfi  =  0  , 

v 

(3) 

2)  cotYxfl  -  -^=-  tanxft  =  0  . 

YP 
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7.  LYAPUNOV-HAUPT  SEQUENCE 


We  now  demonstrate  graphically  that  at  the  end-points  of  the 
Brillouin  zone  the  characteristic  values  of  X  can  be  arranged  in  a 
Lyapunov-Haupt  sequence, [13, 3A],  This  is  an  interesting  and  an  important 
property  of  differential  equations  with  periodic  coefficients,  and  in 
this  relatively  simple  case  it  can  be  demonstrated  explicitly.  We  first 
consider  the  case  of  pure  torsional  mode,  corresponding  to  =  0  . 

Let  yp/p  =  r>  Yt/t  =  T  and  tfi  =  x.  Then  for  c  =  1  ,  eqn.  (6.2) 
can  be  written  as 


1)  Tx  =  kn  -  tan  1 f  tanx  , 


2)  Tx  =  kir  -  tan  tanx  , 


and  for  o  =  -1,  eqn.  (6.3)  takes  the  form 


i)  Tx  =  (2k+l)ir/2  -  tan  lT  tanx  , 


2)  Tx  =  (2k-t-l)TT/2  -  tan-1-^-  tanx  , 


where  in  these  equations  k  =  0,1,2,...  . 

First  consider  the  multi-valued  function 


f(x)  *  kn  -  tan  ^tanx  ,  k  =  0,1,2,... 


(1) 


(2) 


(3) 


m 


tk 


for  a  fixed  value  of  parameter  F  .  The  slope  of  the  function  is 


f'(x)  =  -  f/(cos2x  +  f2sin2x)  , 


(4) 


and  its  curvature  is 


f  =  r(p2 


.  _ sin2x 

(cos2x  +  r2sin2x)z 


(5) 


Since  r  >  0,  the  slope  is  always  negative  and  reaches  a  value  -T  for 

x  =  0 ,ir,2ir, . . . ,  and  a  value  of  -1/F  for  x  =  tt/2,3tt/2,57t/2 . 

For  F  =  1  ,  the  curvature  is  zero.  For  r  ^  1,  the  curvature  is  zero 
when  x  =  0  ,tt/2 ,tt  ,3n/2 , . . .  .  When  T  =  1,  f  (x)  =  kiT-x,  and  represents 
a  straight  line  with  slope  -1  and  zero  curvature.  Further  f (x)  =  kn 

for  x  =  0,71,211,...  and  for  x  =  tt/2 ,3tt/2,5tt/2 .  f(x)  =  (k-  y)ir  . 

For  other  values  of  x  ,  the  function  can  be  readily  computed  and  the 
spectral  lines  are  shown  in  Fig.  (3),  marked  1  .  We  can  similarly 
analyze  the  function 

f  (x)  =  kir  -  tan  tanx, 


and  the  spectral  lines  are  shown  in  Fig.  (3),  marked  2.  The  points  of 
intersection  of  the  set  of  these  two  family  of  curves  with  the 
straight  line  with  slope  T  locates  the  zeros  of  the  two  characteristic 
equations  (6.2)  .  It  is  obvious  from  the  figure  that  the  zeros  of 

these  two  characteristic  equations,  which  correspond  to  periodic  solutions 
with  period  d  ,  form  a  sub-sequence 
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o  =  a  <  n  $  ft-  <  ft0  s  ft,  <  q,  (  ft,  <  . , 

O  ±  L  1  4  3  D 


(6) 


We  can  similarly  plot  the  two  functions  f(x)  corresponding  to 

the  case  o  =  -1  .  In  Fig.  ( 3  ),  these  curves  are  marked  i  and  2  . 

The  intersection  of  these  curves  with  the  straight  line  with  slope  T  , 

locates  the  zeros  of  the  two  characteristic  equations  (6.3)  .  If 

1 » 2 
A 

ft  represents  an  element  of  this  sequence,  then  it  is  obvious  from  the 
figure  that  the  elements  can  be  ordered  in  a  sub-sequence 


ft,  S  ft-  <  ft0  $  ft,  <  ft,  <  ft,  < 

1  2  J  4  5  6 


(7) 


Combining  the  two  sub-sequences,  we  have  the  Lyapunov-Haupt  sequence 


o  *  a  <  fL  $  ft-  <  ft,  $  ft„  <  ft.  s  ft,  <  ft-  s  ft,  < 

012  1  2  J  4  J  4 


(8) 


It  may  be  remarked  that  such  a  sequence  is  a  characteristic  of  a  wide 
class  of  problems  governed  by  differential  -quations  with  periodic 
coefficients. 

We  now  consider  the  general  case  for  the  radial  eigenvalue  k2  >  0. 
In  this  case  the  characteristic  eigenvalues  at  the  end-points  of  the 
Brillouin  zone,  for  o  =  +1  and  o  =  -1  are  given  by 


1)  t£  =  kTT  -  tan  1  tamC  , 


2)  t£ 


kw  -  tan  1  w  tanxC  ; 


(9) 


1)  t X  -  (2k+l)y  -  tan  1  tanrC  , 

2)  t 5  =  (2k+l)y  -  tan  1  tanx^  • 
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Consider  the  first  of  the  four  equations  (9)^»  In  the  notation 
defined  earlier,  this  equation  can  be  rewritten  as 


TYx^Ckt/y)2  =  k/rr  +  f^(x)  , 


(10) 


where  the  function  f^(x)  is  defined  as 


f1(x)  =  -tan  1 | F 


Jix  t)^)2  tanA2_(KT)2 


,  x  >  kt/y 


=  0  , 


x  =  kt/y 


-i  tanh 


-1 


—  tan/x2- (kt ) ^  ,  kt  <  x  <  kt/y 


0  ,  x  =  / (kt) 2+  tt^ 


=  0  ,  x  =  / (kt) 2+  4tt^ 


-i  tanh 


-1 


Tkt/ (1/y)^-1 


,  x  =  kt  ,  Y  <  1 


=  -i  tanh 


-1 


r  Akt/ Y)^-x? 
y  (kt)2-x2 


tanlV (kt)2-x^  ,  0  <  x  < 


=  -i  tanh 


-1 


—  tanhKT 

(Y  ) 


]  ,  x  0  . 


For  t  =  b/a  =  1/3  ,  the  first  radial  eigenvalue  is  k^  =  1.18920.  For 


the  purpose  of  numerical  computations  we  select  as  before 


p/p  *  40,  p/p  =  5/2,  i/a  =  3/5,  I/a  =  3  , 
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which  leads  to 

Y  =  1/4,  r  =  10,  T  =  9 it/20,  t  =  9tt/4,  T  =  5/4  , 

as  the  values  of  the  parameters.  We  now  plot  the  function  kit  +  f^(x)  , 

for  k  =  0,1,2 .  The  points  of  intersection  of  these  curves  with 

the  parabola  Ti/x2-(kt /y) 2  ,  locate  the  zero  of  the  characteristic  equa¬ 
tion  (9)^*  In  Fig.  (4),  the  parabola  and  the  spectral  lines  are  drawn 
for  x  i  kt/y  ,  and  the  spectral  lines  corresponding  to  the  character¬ 
istic  equation  (9)^  are  marked  1  .  The  spectral  lines  for  the  second 
characteristic  equation  are  also  shown  in  the  same  figure,  and  marked  2. 
It  is  obvious  from  the  figure  that  the  zeros  of  the  two  characteristic 
equations,  when  arranged  in  an  increasing  order,  form  a  subsequence, 
which  is  similar  to  the  case  kq  =  0  .  The  second  system  of  character¬ 
istic  equations  are  plotted  similarly  and  are  shown  in  the  same  figure. 
These  spectral  lines  are  marked  £  and  2  ,  respectively.  Again  it  is 
abvious  from  the  figure  that  corresponding  to  the  zeros  ,  there  again 
exists  a  sub-sequence,  similar  to  the  case  =  0  .  Combining  the  two 
subsequences,  we  find  that  the  zeros  of  these  four  equations,  when 


arranged  in  an  ascending  order,  form  the  Lyapunov-Haupt  sequence. 


8.  CO-EXISTENCE  OF  PERIODIC  SOLUTIONS 


To  study  the  problem  of  co-existence  of  periodic  solutions,  we 
first  consider  the  case  of  pure  torsional  mode  when  k  =  0  .  For  a  =  1, 
the  two  characteristic  equations  are 


T  sinx  cosTx  +  cosx  sinTx  =0  , 

T  cosx  sinTx  +  sinx  cosTx  =0 


(1) 


Consider  first  the  case  when  F  =  1  .  In  Fig.  (  3 ) ,  the  spectral  lines 
corresponding  to  F  =  1  are  straight  lines  with  slope  -1.  The  points 
of  intersection  of  these  lines  with  a  line  of  slope  T  are  all  roots  of 
multiplicity  2.  This  also  follows  from  the  fact  that  when  r  =  1  ,  these 
two  equations  are  equivalent  to  one  equation  cos2(l+T)x  =  1  ,  whose 
zeros  are  given  by  =  mr/(l+T)  ,  n  =  0,1,2,.,.  .  These  are  zeros 
of  multiplicity  2,  because  the  derivative  sin2(l+T)x^  =  0  .  Thus, 
when  o=l  and  T  =  1  ,  we  arrive  at  the  ft-subsequence 

o  =  n  <  n  =  <  u,  ■=  n.  <  a  =  ^  <  ...  . 

o  1  2  3  4  5  6 

Similarly  when  o  =  -1  ,  the  two  characteristic  equations  are 

F  sinx  sinTx  -  cosx  cosTx  =  0  , 


T  cosx  cosTx  -  sinx  sinTx  =  0 


(2) 
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When  T  =  1  ,  these  two  equations  are  equivalent  to  one  equation 

cos2(l+T)x  =  -1  ,  whose  zeros  are  x  =  (n+1/2) tt / (1+T)  ,  n  *=  0,1,2,. 

n 

Again  these  are  zeros  of  multiplicity  2,  because  the  derivative 

sin2(l+T)x  =  0  .  Thus  in  this  case  we  arrive  at  the  fi-subsequence 
n 


£2,  =  fL  <  £L  =  n.  <  SL  =  ft,  <  ...  . 

1  2  3  H  j  6 


Combining  the  two,  we  have  in  this  case  the  Lyapunov-Haupt  sequence 


0  —  £2  <  £2,  “  £2*  <  £2,  —  £2«  ^  £2~  =  £2,  <  £2~  ~  £2#  ^  ...  ,  (3) 

ol2123434 


for  T  =  1  . 

We  now  consider  the  case  when  T  ^  1  ,  a  -  1  .  To  find  the  zeros 
of  multiplicity  2,  we  consider  the  product  of  eqns.  (1)^  and  (1)2*  which 
can  be  written  in  the  form 


/F - -  sin2x  sin2Tx  -  2 [cos2 (1+T)x-1]  =  0  . 

/F 


The  first  derivative  of  this  equation  is 


/F  -  —  [cos2x  sin2Tx  +  Tsin2x  cos2Tx]  +  2(1+T)sin2(l+T)x  =  0. 

/F 

(5) 

Evidently,  these  two  equations  will  be  satisfied  if  x  and  T  are  so 
chosen  that 


sin2x  *  0  ,  sin2Tx  =  0  ,  and  cos2(l+T)x  *  1  . 
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The  roots  x  so  determined,  are  roots  of  multiplicity  2  and  are  given 
n 


by 


xn  =  q  V2  ,  T  =  p/q 

(q+p)  =  2n,  q,p,n  =  0,1,2,3,... 


(7) 


We  thus  see  that  there  exist  multiple  roots  only  when  T  is  a  rational 
number.  For  rational  values  of  T  ,  a  line  with  this  slope  will  pass 
through  the  intersection  points  of  the  spectral  lines  corresponding  to 
the  two  equations  (1).  The  coordinates  of  the  intersection  points  are 
(q  y  ,  p  y)  ,  where  at  each  intersection  point  (q+p)  =  2n  .  These 
intersection  points  are  clearly  shown  in  Fig.  (3),  and  all  those 
intersection  points  which  lie  on  the  line  passing  through  the  origin 
with  slope  T  =  p/q  ,  are  zeros  of  multiplicity  2.  Thus  if  (q  it/2, 
p  ir/2)  where  (q+p)  =  2n  be  the  coordinates  of  the  first  zero  of  mul¬ 
tiplicity  2  on  the  line  T  =  p/q  ,  then  the  coordinates  of  higher  zeros 
on  this  line  are  m(q  it/2,  p  tt/2)  ,  where  m(q+p)  =  2n.  As  a  typical 
case,  the  SI- sub  sequence  for  T  =  3/1  is  given  by 


0  =  U  <  SI,  <  <  S2„  =  ft.  <  <  SJ-  *=  0o  <  <  . . .  (8) 

0123456789 

We  now  consider  the  second  case  when  T  ^  1  and  o  =  -1  .  It 
can  similarly  be  shown  that  corresponding  to  equation  (2),  there  exist 
double  roots  if  x  and  T  are  so  chosen  that  they  satisfy  the  equations 


sin2x  =  0  , 


sin2Tx  =  0 


and  cos2(l+T)x  =  -1 


(10) 
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In  this  case,  we  have  roots  of  multiplicity  2  when 

xn  =  q  n/2  ,  T  =  p/q  (11) 

(q+p)  =  2n+l,  q,p,n  =  0, 1*2 ,3, . . . 

The  coordinates  of  the  intersection  point  of  the  family  of  curves  corres¬ 
ponding  to  equation  (2)  are  (q  tt/2,  p  tt/2)  where  (q+p)  =  2n+l  .  Those 
intersection  points  which  lie  on  the  line  with  slope  T  =  p/q  ,  have 
the  coordinates  m(q  tt/2,  p  u/2)  ,  where  m(q+p)  =  2n+l  .  As  a  typical 
case,  the  fi-subsequence  for  T  =  2/1  is 

^  <  Qy  <  =  £210  <  ...  (12) 

Combining  the  two  results,  the  Lyapunov-Haupt  sequence  for  a  typical 
value  of  T  =  p/q  can  now  be  easily  constructed.  Choosing  T  =  3/1  , 
the  sequence  is 

o  =  n  <  n,  <  fL  <  ft,  <  <  £L  <  n,  <  ft,  =  n. 

012123434 

<  <  &7  <  fig  <  n7  =  <  ...  (13) 

Now  consider  the  intersection  of  the  line  with  slope  T  =  p/q  , 
with  the  spectral  lines  corresponding  to  o  =  ±1  .  At  each  intersection 
point  there  exists  at  least  one  periodic  solution,  of  period  d  or  2d  , 
depending  upon  a  m  +1,  or  -1  ,  respectively.  When  the  line  with  slope 
T  -  p/q  intersects  a  lattice-point  with  coordinates  (q  n/2,  p  n/2)  , 
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two  periodic  solutions,  both  of  period  d  or  of  period  2d  ,  co-exist. 
This  is  called  co-existence, implying  existence  of  2  periodic  solutions 
of  the  same  period,  or  coincidence,  implying  existence  of  eigenvalues 
of  double  multiplicity.  Naturally,  coincidence  of  frequencies  implies 
co-existence  of  periodic  solutions,  and  conversely. 

We  now  consider  the  problem  of  co-existence  of  periodic  solutions 
when  the  radial  eigenvalue  k2  >  0  .  The  theory  tells  us  that  when 
o  =  ±1  ,  there  exists,  in  general,  one  periodic  solution  corresponding 
to  every  simple  eigenvalue  of  the  characteristic  equations  (5. 3, 5. 4).  The 
other  linearly  independent  solution  corresponding  to  the  same  eigenvalue 
is,  in  general,  aperiodic.  The  problem  of  co-existence  of  periodic 
solutions  requires  the  existence  of  non-simple  eigenvalues  of  double 
multiplicity.  We  therefore  look  for  double  roots  of  the  characteris¬ 
tic  equation  when  k2  >  0  .  When  o  =  1  ,  k2  >  0  ,  the  characteristic 
equation  (5.3)  has  the  equivalent  form 


SR  -  — 

SK> 


sin2/x2-(Kt)2  sin2T/x2-(KT/y) ^ 

-  2  cos2 (/x2-(kt) 2  +t/x2- (kt/y) 2J  -lj  =  0  ,  (14) 


where 


A  =  IV(x2-(kt/y)2) / (x^-(kt)2)  , 

r  =  Yh/w  »  T  =  yt/t  ,  x  h  t!1 


(15) 


We  consider  the  case  when  1/y  >  1  .  Then  equation  (14)  and  its  first 
derivative  will  simultaneously  be  satisfied  if  x  and  T  are  so  chosen 
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that  they  satisfy  the  equations 

sin2/x2-(Ki)'’  =  0  ,  s  in2lVx2- (kt /y)  *  =  U  , 

(16) 

cos2 (/x^- (kt) 2  +  T/x2-  (kt/y)  2 )  *  1 

The  roots  of  multiplicity  2  are  therefore  given  by 

x  =  / ( qTr / 2)^f  (kt ) 2  T  =  I  ■  ■  ■  ?  y  -  ■<  ,  (17) 

p  +  q  =  2n  ,  q,p,n  =  0,1,2,3, . . . 

We  thus  see  that  there  exist  roots  x^  of  multiplicity  2  if,  in  general, 

T  is  an  appropriately  chosen  irrational  number  satisfying  condition 
( 17 ) 2 •  If  for  a  given  k  ,  there  exist  an  integer  Z  and  parameters 
y  and  t  such  that  (l/y2-l)  (2t/tt)2  =  (q2-i2)/nc2  ,  then  T  is  a 
rational  periodic  decimal  number.  If  we  assume  that  T  is  a  suitably 
chosen  irrational  number,  then  a  parabolic  curve  Tv/x^-(k.t/y)‘  will 
pass  through  a  lattice-point  (x^,  p-n /2 )  of  the  intersecting  family  of 
spectral  lines  corresponding  to  the  two  characteristic  functions  (7.9)^  ^ • 
This  lattice-point  (x^,  pir/2),  p  +  q  =  2n,  is  then  a  zero  of  multiplicity 

2.  Thus  for  y  =  1/4,  t  =  9w/20,  <1  =  1.18920,  x$  =  2.536  633”  will 

be  a  zero  of  multiplicity  2  if  we  choose  t  in  such  a  way  that 

T  ■  p/2.795  981,  where  5  +  p  “  2n.  Similarly  =  3.047  355”  will  be 

a  zero  of  multiplicity  2  for  T  ■  p/4.337916,  where  6  +  p  *  2n.  As  an 
example,  let  q  *  6  and  p  *»  2  so  that  4.337  916T  =  2.  The  parabola 
T/xz-(kt/y)2  now  intersects  the  lattice-point  (xfa,  u)  ,  and  examination 


of  Fig.  (  4  )  reveals  that  we  have  a  subsequence  which  can  be  arranged  in 
an  ascending  order,  except  that  one  of  the  eigenvalue  is  of  multipli¬ 
city  2.  The  structure  of  this  sequence  is  obvious  from  Fig.  (4). 

We  now  consider  the  case  when  o  =  -1  ,  k2  >  0  .  It  can  similarly 
be  shown  that  corresponding  to  the  characteristic  equation  (5.4),  there 
exist  double  roots  if  x  and  T  are  so  chosen  that  they  satisfy  the 
equations 


(18) 


In  this  case  the  roots  of  multiplicity  2  are  given  by 


x  =  /(qu/2)2+  (kt)2  ,  T  =  — 7- 

q  /2 


p  +  q  =  2n  +  1  ,  p,q,n  =  0,1,2,3,... 


The  discussion  concerning  these  roots  of  multiplicity  2,  follows  along 
the  same  lines  as  in  the  previous  case  when  o  =  +1,  and  therefore  hardly 
needs  reiteration.  However,  we  may  repeat  that  existence  of  the  double 
roots  implies  co-existence  of  two  periodic  solutions,  both  of  the  same 
period,  d  or  2d  ,  whichever  the  case  may  be. 
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9.  STABILITY 

To  study  regions  of  lability  and  stability  we  now  plot  graphs  show¬ 
ing  x  versus  T  ,  for  =  0  and  =  1.18920  ,  which  are  shown  in 
Figs.  (5)  and  (6)  ,  respectively. 

Solutions  are  called  stable  and  bounded  when  |h|  «•  2  .  For  values 

of  x  and  T  in  the  hatched  regions  in  these  figures,  all  solutions  are 

stable,  and  the  hatched  regions  are  regions  of  stability.  L'nhatched 

regions  are  regions  of  lability  where  |h|  >2  ,  and  the  solutions  are 

unbounded.  On  the  boundaries  of  hatched  and  unhatched  regions,  |h|  =  2, 

and  therefore  for  a  typical  point  (x,T)  on  the  boundary,  there  exists 

at  least  one  periodic  solution  of  period  d  ,  or  2d  .  At  a  point  of 

intersection,  two  such  periodic  solutions  co-exist.  These  points  of 

intersection,  correspond  to  similar  intersection  points  in  Figs.  (  3  )  and 

(A).  We  note  from  the  stability  diagrams  that  there  is  one  value  of 

T  for  which  Sl^  =  SI^  *  three  values  of  T  for  which  Si^  =  SI^  ,  5  values 

of  T  for  which  SI,.  *  SI^  ,  etc...  ,  no  value  of  T  for  which  fl^  =  , 

2  values  of  T  for  which  SI-  •  Q,  ,4  values  of  T  for  which  SI.  =  SI,  , 

3  4’  50 


etcetera. 
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Fig.  (1)  Floquet's  spectrum  on  an  extended  zone  scheme  for  pure 
torsional  modes  kq  =  0.  Heal  segments  are  shown  as 
full  lines,  imaginary  segments  are  shown  as  dotted 
lines.  Hatched  regions  represent  stopping  bands, 
i  [  >•  i :  (0  -  tt/2(t+t)). 
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Fig.  (2)  Floquet's  spectrum  on  an  extended  zone  scheme  for 

torsional  modes  with  one  radial  node,  =  1.1892. 

Real  segments  are  shown  as  full  lines,  imaginary 
segments  as  dotted  lines.  Hatched  regions  repre¬ 
sent  stopping  bands.  Wave  number  £  =  0  at  point  A, 
and  at  point  B,  7,  *  0.  Brillouin  zone  1  X ] :  (O-tt/2  (t+t)  ) . 


1  5t/(t+t) 


^  =1.18920 


Fig.  C4)  Plot  of  the  multi-valued  functions  f(x)  versus  x  for 
K  =  1.1892,  and  for  x  5  kt/y  =  2.140b6*.  A  typical 
parabolic  curve  t/x?-(kt/7)T  is  also  shown  intersecting 
the  multivalued  function  f(x)  . 


T  for  pure 


Stability  chart  x  versus 
torsional  .odes,  -  "• 

stabliitv  are  shorn  batched  .here  the  sold- 


tions  are  bounded, 
co-exist  at  points 


Periodic  solutions 


ot  intersection. 


*  =  1.18920 


0.5  1 

fig.  (6)  Stability  chart  x  versus  T  for 

torsional  modes  with  one  radial  node, 
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